Spin dynamics in the regime of hopping conductivity 
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We consider spin dynamics in the impurity band of a semiconductor with spin-split spectrum. Due 
to the splitting, phonon-assisted hops from one impurity to another are accompanied by rotation of 
the electron spin, which leads to spin relaxation. The system is strongly inhomogeneous because of 
exponential variation of hopping times. However, at very small couplings an electron diffuses over a 
distance exceeding the characteristic scale of the inhomogeneity during the time of spin relaxation, 
so one can introduce an averaged spin relaxation rate. At larger values of coupling the system 
is effectively divided into two subsystems: the one where relaxation is very fast and another one 
where relaxation is rather slow. In this case, spin decays due to escape of the electrons from one 
subsystem to another. As a result, the spin dynamics is non-exponential and hardly depends on 
spin-orbit coupling. 

PACS numbers: 71.70.Ej, 73.63.-b, 75.40.Gb, 85.75.-d 



Spin dynamics in semiconductors has attracted much 
attention in the last decades (for review, see Ref. [l|). A 
great number of publications were devoted to the study 
of the spin relaxation in the metallic regime, where it 
is usually dominated by the Dyakonov-Perel mechanism. 
According to this mechanism, the spin relaxation rate is 
proportional to the diffusion coefficient [2] . Recently, the 
idea was put forward [3j that this proportionality goes 
beyond the metallic regime and stays valid in the regime 
of hopping conductivity: 
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Here r^ 1 ' is the spin relaxation time, L$ ~ v/fl, v is 
the tunnelling velocity, and Q is the typical value of spin 
precession frequency for an electron moving under the 
barrier. 

In this paper we discuss spin dynamics in the hopping 
regime and show that a number of spin relaxation mech- 
anisms are realized in addition to the one discussed in 
Ref. 0]. The hopping regime takes place at small tem- 
peratures T <C W (here W is the width of the impu- 
rity band). In the opposite limit T ^> W the transport 
is governed by activation to the mobility edge. How- 
ever, methodologically it is convenient to start from the 
modelling case when T ^> W but transitions to the con- 
duction band are forbidden. This model is rather simple 
because the probability of a hop does not depend on tem- 
perature. On the other hand, it captures basic physics of 
the problem and can be easily generalized to the real sit- 
uation. We neglect electron-electron interaction and as- 
sume that the distance between impurities is much larger 
than the Bohr radius a: na d <C 1 (here n is the impurity 
concentration). We also assume that Lg ^S> a which is 
usually the case. 

Due to the spin-orbit coupling, phonon-assisted hops 
are accompanied by rotation of the electron spin. The ro- 
tation angle can be written as tp = (p +Scp, where tp cor- 
responds to the quasiclassical approximation, which im- 



plies that electron moves under barrier along the straight 
line connecting the impurities and 5<p is the correction 
which is due to quantum uncertainty of the trajectory 
(this correction cannot be neglected when we consider 
non-typical hops over the distances of the order of sev- 
eral a). The angle ip is given by ipo = 0(v)t, where 
fl(v) is proportional to spin-orbit spectrum splitting, 
v = nfi/raa is the tunnelling velocity, t — r/v, a is the 
Bohr radius, m is the electron mass, r is the hopping 
length, and n is the unit vector in the direction of the 
hop. In the 3D case fl x ~ v x (Vy — v^) and the other com- 
ponents are given by cyclical transmutations [4|. In the 
2D case fi = &v x + hv y , where constant vectors a and 
b depend on the quantum well orientation and the de- 
gree of asymmetry of the confining potential The 
characteristic value of ipo is given by 



ipo ~ r/ L t 



(2) 



Random hops between impurities lead to spin relax- 
ation. Let us assume that at t = the homogeneous 
distribution of the spin density So is created in the sys- 
tem (for example, by means of optical orientation) and 
study the dynamics of the averaged spin density S(i). 
The system under discussion is strongly inhomogeneous 
due to exponential variation of hopping times. Indeed, 
the waiting time Tjj of a hop between impurities i and j 
varies exponentially with the distance between impu- 
rities: = to exp(2Ly/a), where tq is the waiting time 
for the impurities separated by the distance of the order 
of a (here we took into account that T 3> W). We as- 
sume 0| that impurities i and j are effectively connected 
if t^ < t, or, equivalently, 

L i:j < L(t)=a/2ln(t/T ), 

and disconnected if > t (Lij > L(t j). For L(t) larger 
than a certain critical length Lh ~ n^ 1 / d , the connected 
impurities form infinite cluster. The fraction of impuri- 
ties P(£) belonging to the infinite cluster increases with 
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increasing of t as shown in Fig.l. Here 

e = £(t) = L(t)/L h = (a/2L h ) ln(t/r ). (3) 

Near the critical point (L(t) — Lh) the share of impurities 
belonging to the infinite cluster P behaves as a power 
function: P = (£ — I) 73 , (where (3 is critical index), and 
with the further increase of £ it is quickly saturated 0] ■ 
The exponential regime of the spin relaxation S(t) = 
Soexp(— t/rg) with an averaged rate 1/ts is realized at 
sufficiently small spin-orbit coupling, when an electron 
has enough time to visit almost any possible impurity 
configuration during the time t$- This condition is sat- 
isfied when 

1 - P[Z(t s )] « 1. (4) 

There are two contributions to 1/t<j. The first one is due 
to electron diffusion on the scales larger than the scale 
of inhomogeneity. This contribution is given by Eq. ([1]) 
Q. Another contribution, neglected in Ref. [3], comes 
from the spin relaxation on non-typical clusters with the 
distances between impurities on the order of several a. 
The most likely non-typical configurations of impurities 
are pairs and triangles. An electron captured on a non- 
typical complex makes many hops in a small region of 
space before it leaves it. The shift in space due to these 
hops is limited by the size of the complex, while the de- 
viation of its spin increases with the number of hops. 
Therefore, such complexes do not contribute to the dif- 
fusion process but they might give essential contribution 
to the spin relaxation. 

The most probable non-typical configuration is a pair 
of impurities. In quasi-classical approximation the con- 
tribution of pairs to spin relaxation is equal to zero, be- 
cause the rotation of the spin during the hop from one 
impurity to another is compensated by the rotation dur- 
ing the hope in the opposite direction. However, when 
the distance between impurities is of the order of a, the 
deviation of the trajectory from the quasi-classical one 
due to quantum uncertainty is of the order of the size of 
the quasiclassical trajectory itself. This implies that the 
correction Sip to the rotation angle is comparable with ip 
(i.e. Sip ~ ipo ~ a/Lg <C 1). This is indeed the case in 
3D semiconductors. In the 2D case the situation is more 
subtle. In this case the rotation angle in the first order of 
X/Ls depends only on the shift of the electron in space 
ip = J Sldt = ar x + br y (here the integration is taken 
over the electron trajectory and we took into account 
that small rotations commute) and does not depend on 
a particular trajectory. As a result, the correction to 
the main contribution appears only in the next orders in 
1/Ls- Sip ~ (a/Ls) m , where to > 2. One can show that 
?7i = 2 for the components of the spin which are perpen- 
dicular to [a x b] and to — 3 for the parallel component. 
The exact value of this correction depends on the wave 
vector of the phonon that was emitted or absorbed dur- 



ing the hop [8j . Therefore, the rotation angle is a random 
quantity. 

In general case, the angle of the rotation after passing 
the closed trajectory (that includes the hop from the first 
impurity to the second one and the hop in the opposite 
direction) can be written as: 

Sip ~ (a/L s ) m , (5) 

where to = 1 in 3D case and m = 2, 3 in 2D case. 

After rotation by the angle Sip, spin projection on the 
original direction decreases on the average by Sip 2 . So, 
one can expect that a relaxation rate on a single pair is 
given by (a / L s) 2m (1 / Vo) . However, this estimate does 
not take into account the effects of level repulsion. The 
repulsion energy is on the order of AE ~ (e 2 /«;a)e~ r l a . 
A non-typical pair can participate in the spin relax- 
ation only if AE < T. Therefore, the minimal dis- 
tance between impurities in a non-typical pair is given 
by r* > a\n(Ka/Te 2 ). The hop waiting time for such a 
pair is given by t exp(2r*/a) ~ r (e 2 /Tna) 2 , this leads 
to an additional factor (T/ta/e 2 ) 2 <C I in the spin relax- 
ation rate. As a result, the rate of the spin relaxation on 
a single non-typical pair reads 

l/T^(a/L s ) 2m (l/r )(T K a/e 2 ) 2 . (6) 

The contribution of non-typical pairs to the spin relax- 
ation rate is proportional to 1/rg and to the share of 
the time that electron spends on non-typical pairs. This 
share is proportional to the ratio of the concentration of 
non-typical pairs n* to the total impurity concentration: 
n* jn ~ (a/Lh) . Hence, the corresponding relaxation 
rate is given by 

1/rf ~ {a/L h ) d {a/L s ) 2m (l/r Q )(T K a/e 2 ) 2 . (7) 

Both mechanisms described above are additive and the 
spin relaxation rate is: 

l/r s = l/4 1) +l/4 2) . (8) 

Let us compare these two mechanisms. In the case un- 
der discussion, when T ^> W, the diffusion coefficient 
entering Eq. ([!]) is determined by the hops over the crit- 
ical length: D ~ cxp(-2L h /a) [7]. Therefore, 
decreases exponentially with the average distance be- 

(2) 

tween impurities, while I/t^ decreases according to slow 
power law 1/Lf (see Eq. (0). In the 3D case the de- 
pendence on Ls in Eq. is the same as in Eq. (fT|). 
Therefore, in the 3D case, l/4 2) > 1/4 1} (since we as- 

(2) 

sumed that Lh 3> a). In the 2D case, 1/Tg decreases 
with Ls faster than I/t^ 1 '. Thus, at very small spin- 
orbit coupling the main contribution is given by Eq. |T]). 
However, the crossover to the regime of non-typical pairs 
takes place already at exponentially small value of the 
coupling. 
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In deriving Eq. we assumed that during the pro- 
cess of spin relaxation an electron has enough time to 
visit many non-typical pairs which are separated by a 
distance on the order of [n*] -1 / d = L^L^/a) (this is the 
length that plays the role of inhomogcneity scale in this 
regime) . With the increasing of the strength of spin-orbit 
coupling this condition fails because spin relaxation time 
decreases and there appear finite clusters that are effec- 
tively separated from the rest of the system on the time 
scale on the order of ts- Since the concentration of non- 
typical pairs n* is small, most of the isolated clusters do 
not contain them. This means that the spin of electrons 
captured on the isolated clusters does not relax until they 
escape to the infinite cluster where the relaxation is much 
faster (it is still determined by Eq. ([8])). Therefore, the 
infinite cluster can be considered as a "black hole" for 
electron spin. The spin dynamics in this regime is given 
by 



of electrons captured on such " traps" : 



S(t) = S [l-P(t)], 



(9) 



where P(t) = P[L(t)/L h ] = P[(a/2L h )1n(t/r )]. Hence, 
the spin polarization is a function of a single dimension- 
less parameter (a/2Lh) ln(t/ro). It worth noting that the 
dynamics of the spin relaxation in this regime does not 
depend on the strength of spin-orbit coupling. 
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FIG. 1: Fraction of impurities belonging to the infinite cluster 
as a function of the maximal distance L between connected 
impurities. In the vicinity of the critical length Lh, function 
P[L/Lh] increases as [L/Lh — l]' 3 , where is a critical index. 



With the further increase of the coupling, the spin re- 
laxation on the finite clusters comes into play. Finally, 
at extremely high coupling the spin rotates by a large 
angle after one typical hop. In this case, the spin relaxes 
on any finite clusters except isolated impurities separated 
from the rest of the system by distances larger than L(t). 
The spin at the moment t is proportional to the number 



S(t) = Sq exp 



J->h T 



(10) 



where d = 2, 3 is dimensionality of the system and Cd ~ 1 
is a numerical coefficient. It is seen from this equation 
that just as in the case of intermediate coupling the spin 
dynamics does not depend on the spin-orbit strength and 
determined by the only parameter (a/Lh) ln(i/ro). 

Next, we discuss qualitatively the generalization of 
the approach discussed above to the case of small tem- 
peratures. In this case the waiting time is given by 
T ij = T o CX P (Li j I a + Eij/T), where E^ is the energy 
difference between states i and j. The minimization of 
this exponent yields the optimal values of hopping length 
Lij ~ Lh(T) and hopping energy £Jy ~ Eh(T) which are 
given 



by 



E h (T) ^ L h (T) 
T ~ a 



(9Ta d ) 



d\-i/(d+i) 



(11) 



(these equations are applicable when T < Wa/Lh). 
The optimal length Lh(T) is of the order of the av- 
erage distance between impurities lying in the energy 
band of width E^T) (optimal energy band). The states 
outside the band of width Eh(T) near the Fermi level 
are not accessible to polarized electrons. As a conse- 
quence, the equations derived above are still valid with 
two modifications. First, one should make a replacement 
Lh —> Lh(T). Second, the effective dimension of the sys- 
tem should be increased d — > d + 1, which takes into 
account the motion of the electron along energy axis. 

To conclude, the theory of the spin relaxation of the 
electrons in the impurity band was developed and a num- 
ber of different relaxation regimes were predicted. At 
weak spin-orbit couplings spin relaxation is exponential 
with the rate of spin relaxation given by two contribu- 
tions: the first one comes from the relaxation on non- 
typical pairs of impurities and the second one is due to 
electron diffusion on the scales larger than the scale of in- 
homogeneity of the system. At stronger couplings, spin 
relaxation is due to escapes of the electrons from finite 
clusters to the infinite one. In this case, the law of spin 
relaxation does not depend on the spin-orbit coupling. 
At very large couplings the spin relaxation is due to es- 
cape from the spin traps. 
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